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We analyze low-energy scattering for arbitrary short-range interactions plus an 
attractive 1/r 6 tail. We derive the constraints of causality and unitarity and find 
that the van der Waals length scale dominates over parameters characterizing the 
short-distance physics of the interaction. This separation of scales suggests a sepa- 
rate universality class for physics characterizing interactions with an attractive 1/r 6 
tail. We argue that a similar universality class exists for any attractive potential 
l/r a for a > 2. We also discuss the extension to multichannel systems near a mag- 
netic Feshbach resonance. We discuss the implications for effective field theory with 
attractive singular power-law tails. 



I. INTRODUCTION 



Low-energy universality appears when there is a large separation between the short- 
distance scale of the interaction and the physically relevant long-distance scales. Some 
well-known examples include the unitarity limit of two-component fermions [HE] and the 
Efimov effect in three-body and four-body systems [oTfTo]. See Refs. JTUlIITj for reviews of the 
subject and literature. There have been many theoretical studies of low-energy phenomena 
and universality for interactions with finite range. These studies have direct applications 
to nuclear physics systems such as cold dilute neutron matter or light nuclei such as the 
triton and alpha particle. To a good approximation, the van der Waals interactions between 
alkali-metal atoms can also be treated as a finite-range interaction. 

However, there are some differences. For potentials with an attractive l/r a tail and a > 2, 
the s-wave scattering phase shift near threshold has been formulated in Ref. [18j. For a > 3, 
the modified scattering parameters for an s-wave Feshbach resonance were derived in Ref. [TH] 
using coupled-channel calculations. Analytical expressions for the s-wave scattering length 
and effective range for two neutral atoms and a = 6 have been derived in Ref. [20] • However, 
the applicability of the effective range theory is limited for interactions with attractive tails. 
In order to define the scattering length for angular momentum L > 2 and the effective range 
for L > 1, a modified version of effective range theory known as quantum-defect theory is 
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needed [211 [22]. Furthermore, scattering parameters of magnetically tunable multichannel 
systems have been studied in the context of multichannel quantum-defect theory [231 121] • 
See Ref. |25J for a very recent development of multichannel quantum-defect theory for higher 
partial waves. There is also growing empirical evidence that there exists a new type of low- 
energy universality that ties together all interactions with an attractive 1/r 6 tail. This might 
seem surprising since there is no such analogous behavior for interactions with a Coulomb 
tail. In this paper we derive the theoretical foundations for this van der Waals universality 
at low energies by studying the near-threshold behavior and the constraints of causality. 
We also show that this universality extends to any power-law interaction l/r a with a > 2 
in any number of dimensions. Our analysis applies to energy-independent interactions. We 
first consider a single scattering channel but then also consider multichannel systems near a 
magnetic Feshbach resonance. The full analysis for the multichannel problem off resonance 
will be discussed in future publications. 

In our analysis we assume that the two-body potential has a long-distance attractive tail 
of the form —C^/r 6 . We define the van der Waals length scale, (3$, as 

Aj = (2/xC 6 )*, (1.1) 

where /i is the reduced mass of the scattering particles. For simplicity we use atomic units 
(a.u.) throughout our discussion. So, in particular, we set h = 1. In Refs. [2SJ [2Z] it was 
noticed that an approximate universal relationship exists between the effective range and 
inverse scattering length for s-wave scattering in many different pairs of scattering alkali- 
metal atoms. If we write Aq as the scattering length and Rq as the effective range, the 
relation is 

/ Wl/4) 2 _ Aft 8nPl 
3tt 3A 3r(l/4) 2 A2 

This approximate relation becomes exact for a pure —Cq/t 6 potential. What is surprising 
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about Eq. (1.2) is that the van der Waals length (3q dominates over other length scales which 
characterize the short- distance repulsive force between alkali-metal atoms. This approximate 
universality suggests there is some separation of scales between the van der Waals length /3q 
and the length scales of the short-range forces. This separation of scales will become more 
transparent later in our analysis when we determine the coefficients of the short-range K 
matrix. It would be useful to exploit the separation of scales as an effective field theory 
with an explicit van der Waals tail plus contact interactions. In this paper, we discuss the 
constraints on such a van der Waals effective field theory. 

We note that a similar dominance of the van der Waals length /3 6 has been discovered 
for the three-body parameter in the Efimov effect [2"8H3"U] . In this paper we focus only on 
two-body systems. However, our analysis should be useful in developing the foundations for 
van der Waals effective field theory. This in turn could be used to investigate the Efimov 
effect and other low-energy phenomena in a model-independent way. An extension of our 
analysis may be useful to understand the recently observed universality of the three-body 
parameter for narrow Feshbach resonances [31j . 

The organization of our paper is as follows. We first discuss the connection between 
causality bounds and effective field theory. Next we consider asymptotic solutions of the 
Schrodinger equation. After that, we derive causality bounds for the short-range K matrix 
and consider the impact of these results on van der Waals effective field theory. Then 
we discuss quantum-defect theory and calculate causal ranges for several examples of single- 
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channel s-wave scattering in alkali-metal atoms. We also consider the constraints of causality 
near magnetic Feshbach resonances. We then conclude with a summary and discussion. 



II. CAUSALITY BOUNDS AND EFFECTIVE FIELD THEORY 

For an effective field theory with local contact interactions, the range of the interactions is 
controlled by the momentum cutoff scale. Problems with convergence can occur if the cutoff 
scale is set higher than the scale of the new physics not described by the effective theory. It 
is useful to have a quantitative measure of when problems may or may not appear, and this 
is where the causality bound provides a useful diagnostic tool. For each scattering channel 
we use the physical scattering parameters to compute a quantity called the causal range, 
R b , which is the minimum range for the interactions consistent with the requirements of 
causality and unitarity. For any fixed cutoff scale, the causality bound marks a branch cut 
of the effective theory when viewed as a function of physical scattering parameters [32] 133] . 
The coupling constants of the effective theory become complex for scattering parameters 
violating the causality bound. These branch cuts do not appear in perturbation theory, but 
they can spoil the convergence pattern of the perturbative expansion. 

Wigner was the first to recognize the constraints of causality and unitarity for two-body 
scattering with finite-range interactions [31]. The time delay of a scattered wave packet is 
given by the energy derivative of the phase shift, 

* = 2*. (2.1) 

It is clear that the incoming wave packet must first reach the interacting region before the 
outgoing wave packet can leave. So the causality bound can be viewed as a lower bound 
on the time delay, At. When applied to wave packets near threshold, the causality bound 
becomes an upper bound on the effective range parameter. Phillips and Cohen derived 
this bound for s-wave scattering with finite-range interactions [35J. Some constraints on 
nucleon-nucleon scattering and the chiral two-pion exchange potential were considered in 
Ref. [36J, and relations between the scattering length and effective range have been explored 
for one-boson exchange potentials [37]. As mentioned above, the same authors studied 
the relationship between the scattering length and effective range for the van der Waals 
interaction [2S1 EZ] 

In Refs. [3B1 |3H] the causality and unitarity bounds for finite-range interactions were 
extended to an arbitrary number of space-time dimensions or value of angular momentum. 
A complementary discussion based upon conformal symmetry and scaling dimensions can 
be found in Ref. [10]. Coupled-channel systems with partial- wave mixing were first studied 
in Ref. |32j, and the interactions with attractive and repulsive Coulomb tails were first 
considered in Ref. 1 33 1 . 



III. ASYMPTOTIC SOLUTIONS OF THE SCHRODINGER EQUATION 

We consider a system of two spinless particles interacting via a spherically symmetric 
potential in the center-of-mass frame. As noted in the Introduction, we use atomic units 
where h = 1. The total wave function in the relative coordinate r can be separated into the 
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radial and the angular parts as 



L,M\ 



R^(r)Y L M (f), 



(3.1) 



where is the magnitude of the spatial momentum and (j) are the spherical harmonics 
We define the rescaled radial wave function l/i k Hr) as 



U<?\r) 



rRf{r) 



(3.2) 



We first review the case where the only interaction between the two particles has a finite 
range, R. This means that the interactions are exactly zero when the two particles exceed 
a distance R. Let \x denote the reduced mass of the two-body scattering system. The radial 
Schrodinger equation for a scattering state with energy E = k 2 /(2jj) is then 



d 2 l(l + r 



dr 2 



+ k 2 



Ui k \r) 



R 



2/x / dr'W(r,r')U^\r') 



(3.3) 



We have written the interaction as a rotationally invariant operator with real kernel W(r, r') 
to avoid any assumption regarding the locality or nonlocality of the potential. In our anal- 
ysis we consider only interactions which are energy independent. The finite-range condition 
implies that W(r,r') = if r > R or r' > R. We assume that the interaction is sufficiently 
well behaved at the origin to admit a regular solution. This assumption imposes the re- 
striction that at short distances the potential is not too singular such that the radial wave 
function satisfies the regularity condition 



lim u\ k \r)^U) K> (r) 
L v ' dr L v ' 



d 



r(*)/ 



0. 



(3.4) 



In Ref. 

that 



it is proven that this condition is fulfilled by a class of potentials V(r) provided 



f R , 

/ r \V(r )| dr' < oo . 
Jo 



(3.5) 



We choose a normalization such that, for r > R, the radial wave function has the form 



Ui k \r) 



7rr 




cot SL(k)J L+ i(kr) 



N L+ i 2 (kr) 



(3.6) 



where J L+ i(kr) and N L+ i(kr) are the Bessel functions of the first and second kind, and 
<5l(A;) is the scattering phase shift. 

We now discuss the main problem of interest where the interactions have a long-range 
van der Waals tail. In addition to the non-singular finite-range interaction parameterized 
by W(r, r'), we assume that there is a long-range local potential —Cq/t 6 for r > R. The 
van der Waals length scale fi§ was defined in Eq. (1.1). The radial Schrodinger equation is 



dr A r z r° 



R 



2fi / dr'W(r,r')U {k \r') 



(3.7) 



The step function 6{r — R) cuts off the long-range potential at distances less than R. This 



ensures that we satisfy the regularity condition in Eq. (3.4) and avoids mathematical prob 



lems associated with unregulated singular potentials |42j . The general form of the solutions 



for Eq. (3.7) has been discussed by Gao in Ref. 
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In order to simplify some of the more lengthy expressions to follow, we introduce dimen- 
sionless rescaled variables r s = r//3 6 , k s = f3 6 k, and p s = l/(2r 2 ). In the outer region, r > R, 
the Schrodinger equation reduces to 



L (L + 1) , PI + k2 



or 



dr 



dr 2 s 



+ 



6 







L{L + 1) 



(3.8) 
(3.9) 

using the formalism 

The van der Waals wave functions Fl and Gl are linearly independent solutions of 



Ul k \r) 



0. 



The exact solutions for Eq. (3.9) have been studied in detail in Ref. 
of quantum-defect theory 



Eq. (3.9). In order to write these out we first need several functions defined in Appendix A 



The van der Waals wave functions Fl and Gl can be written as summations of Bessel 
functions, 



F L (k,r) 



,1/2 



Xl(k s ) + Yl(k s 



m=— oo 



X L {k s ) ^ b rn(ks) Ju+m (ps) 

oo 

- Y L (k s ) b m (k s ) N u+m {p s 



, (3.10) 



G L {k,r) 



Xl(k s ) + Yl(k s 



X L {k s ) ^2 h m{k s )N v+m (p s ) 

m=—oo 

oo 

+ Y L (k) 22 b rn(k s ) J u+m (p s 



m=—oo 



. (3.11) 



The function Xl is defined in Eq. (All), and Yl is defined in Eq. (A12). For m > the 



function b m is given in Eq. (A6), while b_ m is given in Eq. (A7). The offset v appearing 



in the order of the Bessel functions is given by the solution of Eq. (A10) in Appendix A 



For notational convenience, however, we omit writing the explicit k s dependence of v. Let 
us define 6 L short \k) to be the phase shift of the van der Waals wave functions due to the 
scattering from the short-range interaction. The normalization of U L k \r) is chosen so that, 
for r > R, 

U L k \r) = F L {k,r) - tanS L shovt \k) G L {k,r) . (3.12) 

Our van der Waals wave functions are related to the functions /£° and g L ° defined of Ref. [13] 
by the normalization factors F^ = /£ / \/2 and Gl = —g L / V%- Henceforth, we write all 
expressions in terms of the short-range reaction matrix 



Kl = tan 5 



(short) 



which is related to the short-range scattering matrix via 



Si 



e 2iS L 



(short) 



i-K L 
i + K l 



(3.13) 



(3.14) 
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For any finite-range interaction, Kl is analytic in k 2 and can be calculated by matching 
solutions for r < R and r > R at the boundary. It can be written in compact form as 



Wi uf\F^) 



W(U 



(k) r (k), 

L i i 



(3.15) 



r=R 



where Uj/ is the solution of Eq. (3.7) that is regular at the origin, and W denotes the 
Wronskian of two functions, 

W(f,g) = fg'-f'g. 



IV. CAUSALITY BOUNDS FOR SHORT-RANGE ^T-MATRIX K L 

In this section we derive causality bounds for the short-range K matrix K^. For this we 
need to expand the wave function Uf/{r) in powers of k 2 . The steps we follow are analogous 
to those used in Refs. [32J E31 ESI EH] • We first expand K L , 

oo 

K L = tan ort) (*)=]£ K L . 2n k 2n . (4.1) 

n=0 

The first two terms K^q and K^ 2 are analogous to the inverse scattering length and effective 
range parameters in the usual effective range expansion. The higher-order terms can be 
regarded as analogs of the shape parameters. Next we expand the van der Waals wave 
functions in powers of k 2 , 

F L (k, r) = f L>0 (r) + f Lj2 (r) k 2 + 0(k 4 ) , (4.2) 

G L (k, r) = g L>0 {r) + g Lt2 {r) k 2 + 0(k A ) . (4.3) 
In the following, we define 

v = \(2L + l), 

which corresponds to the value of v at threshold. Using the low-energy expansions in Ap- 



pendix Bl we find that the coefficients in Eq. (4.2) are 



h,o(r) = r]l 2 J V0 (p s ) (4.4) 

and 

T(u )T(2u - 1) (3 2 1/2[ - , 

h > 2{r) = r(, + i)r ( 2, ) i6 r - ^ (P.) + (Ps) ] 

r(v )T(2v + 1) /3 6 2 1/2f , 

:V [J^+i (Ps) ~ N U0 (p s ) \ . (4.5) 



T(u + l)r(2z/ + 2) 16 



Similarly, the coefficients in Eq. (4.3) are 

g L>0 (r) = r]/ 2 N U0 (p s ) (4.6) 
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and 



r(i/o)r(2m - i) % 1/2[ u 



Using Eq. (3.12), we can now express Uj/{r) as an expansion in powers of k 2 . For r > R, 
we have 

ut\r) = f Lfi {r)-K L ^ Lfi {r) 

+ k 2 [f L , 2 (r) - K Lfi g L ,2(r) - K Lj2 g Lfi (r)] + 0(k 4 ) . (4.8) 

We now consider two solutions of the Schrodinger equation (r) and U^ b \r) with 
momenta k a and kb, respectively. We have 



uf a \r) = 2p\ dr'W(r,r')U^ a, (r'), (4.9) 



/\Tr(fco)/ /\ 



<ir 2 



L(L + 1) , A 



+ ^fl(r - i?) + /c 2 



f/f b) (r) = 2/x / dr'W^rOf/fV) 



(4.10) 



Following the same steps as in Eq. (31)-(36) of Ref. [39J, we obtain the Wronskian integral 
formula 

W[U { L kb \u { L ka) ](r) 



h 2 — h 2 



U^ a \r')U^ b \r')dr' , 



(4.11) 



for any r > R. Using Eq. (4.8) for momenta k a and k^ we find 

W[U?\u^}(r) 



1.2 U2 



W[f L>2 , f L>0 ](r) - K L>0 {W[g L>2 , f Lfi ](r) + W[f L , 2 , g L>0 }(r)} 



(4.12) 



+ Kl i0 W\g L , 2 , g L , ](r) - K Lt2 W[g Lfi , f L ,o)(r) + 0(k 2 a , k 2 h ) . 

In the Wronskian integral formula ( |4.11 ), we set k a = and take the limit kb — > 0. With 
the wave function at zero energy written as the result is 



K 



L.2 



b L (r) 



71 



dr' . 



where 



TV, 



7T 



b L (r) =-W[f L , 2 J L , ](r) + -Ki fi W[g L , 2 ,g Lfi ](r) 



- jK L , {W[g L , 2 , f L ,o](r) + W[f L>2 , g L , ](r)} . 

The Wronskians appearing in Eq. ( |4.14[ ) can be written out explicitly as 

PeP, 



(4.13) 



(4.14) 



W[f L , 2 ,f L , ](r) 



+ 
+ 



16^o (2u - 1) 
16v (2is + 1) 



4(2^,-1) (2i/ + l) 



[Jw -2(Ps) Ju (ps) ~ Ju -l(ps)] 
[J V0+2 {Ps) Ju (Ps) - Ju 0+ l(ps)] 



Ju +l{Ps) Jvo-l(ps) - Jlips) + 



(4.15) 
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W[g Li2 ,g L ,o\(r) 



P&Ps 



+ 
+ 



16z/ (2z/ - 



16z/ (2z/ + 1) 

Peps 

4(2^o-l) (2i/ + l) 



[A^ 0+2 (p s ) iV^(p s )-< +1 (p s )] 



A^ 0+1 (p s ) N U0 ^{p s ) - N*(p s ) + 



(4.16) 



and 



W[^ ) 2ji,o](r) = W[/ i)2 ,^,o](r) 

= 1 ft TY{ J "0-lO S ) [^0+1 (P*) - ^0-1 (Ps)] - N U0 ( Ps ) [Jy (ps) ~ Ju -2(Ps)] } 

16z/ [2u - 1) 



16^(2^ + 1) 



{ J^ 0+ i(p s ) [^ 0+1 (p s ) - JV^_i(p,)] - N uo {p s ) [J U0+2 (p s ) - ^ (P S )] } • 

(4.17) 



The fact that the integral on the right-hand side of Eq. (|4.13|) is positive semidefinite sets 

(4.18) 



an upper bound on the short-range parameter Kl,2- We find that 

K L , 2 < b L (r) 
for any r > R. 



V. IMPACT ON EFFECTIVE FIELD THEORY 

In this section we discuss the impact of our causality bounds for an effective field theory 
with short-range interactions and an attractive 1/r 6 tail. In Fig. [l] we plot the L = 
Wronskians W[f ,2, /o,o], W [g , 2 , go,o], and W [5-0,2, /o,o] for /3 6 = 50 (a.u.). Figures |2j and [3] 
show the analogous plots for L = 1 and L = 2, respectively. 

We note that all of the Wronskian functions in Figs. [TJ [2j and [3] vanish in the limit r — > 0. 
This stands in clear contrast to what one finds for purely finite-range interactions [38, 39J. 
In that case the effective range parameter, r L , satisfies the upper bound 



r L < bf e {r) 

where the function b { [ ee (r) is [cf. Eq. (60) in Ref. 



1 free / 



2r(L-|)r(L 



4 1 



7T 

f \ 2 



-2L+1 



+ 



L + \a L V2 
2tt 



r(L + |)r(L + 



r \ 2L+3 
2/ 



and is the scattering length. Near r = the behavior of b[ ee (r) is 



Lfree / 



2r(L-|)r(L 



7T 



-2L+1 



+ 0(r 2 



(5.1) 



(5.2) 



(5.3) 
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FIG. 1: (Color online) Plot of W[f 0,2, fo,o](r), W[g ^,go,o}(r), and W[g ,2, fo,o](r) as a function of r for 

L = and /3 6 = 50 (a.u.). 




FIG. 2: (Color online) Plot of W[/i,2, /i,o]W, W[fl'i,2,Si,o]('"), and W[gx,2> fi,o]{r) as a function of r for 

L = 1 and /3e = 50 (a.u.). 



We see that b f [ ee (r) diverges to negative infinity as r — > for L > 1. The causality bound 
on ^ also drives to negative infinity for L > 1, 

r ,<- 2r(L -f L+ ^ GPW). 

For an effective field theory with local contact interactions, the range of the interactions 
are controlled by the momentum cutoff scale. No matter the values for ai and ri, it is not 
possible to take the momentum cutoff scale arbitrarily high without violating the causality 
bound for channels with angular momentum L > 1. For finite-range interactions with an 
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FIG. 3: (Color online) Plot of W[f 2 . 2 , /j.olW, W[g 2 ,2,g2,o](r), and W[g 2 ,2, f2,o)(r) as a function of r for 

L = 2 and f3 6 = 50 (a.u.). 

additional attractive or repulsive Coulomb tail, one finds the same leading behavior [33] 1 

b „ r) = _ 2r(L-l)r(L + i) ( ,p + 

i.e., the only difference in the causality bound relation for the Coulomb-modified effective 
range is the subleading 0(r~ 2L ) pole term which is absent in the purely finite-range case. 
Hence, also for an effective field theory with contact interactions and long-range Coulomb 
tail, it is not possible to take the momentum cutoff scale arbitrarily high for L > 1 without 
violating the causality bound. 

There is no such divergence in &t(r) at r = for the attractive 1/r 6 interaction. For an 
effective field theory with contact interactions and van der Waals tail, the causality bound 
does not impose convergence problems as long as Kl,2 is less than or equal to zero. This 
holds true for any L. There is no constraint from causality and unitarity preventing one 
from taking the cutoff momentum to be arbitrarily large. The key difference between the 
van der Waals interaction and the Coulomb interaction is that, when extended all the way to 
the origin, the attractive 1/r 6 interaction is singular and the spectrum is unbounded below. 
An essential singularity appears at r = 0, and both van der Waals wave functions F L and 
Gl vanish at the origin. 

These exact same features appear in any attractive l/r a interaction for a > 2 in any 
number of spatial dimensions. The same can be said about an attractive 1/r 2 interaction 
when the coupling constant is strong enough to form bound states. The key point is that in 
the zero-range limit of these attractive singular potentials, the spectrum of bound states ex- 
tends to arbitrarily large negative energies. As a consequence, the scattering wave functions 



1 To get this analogy, we use here the normalization of the Coulomb-modified effective range expansion 
found in Eq. (28) of Ref. [33] and insert it in Eqs. (64), (A. 6), and (A. 7) of the same paper, which give the 
explicit expressions for the Coulomb- modified causality bound functions for L = 0,1,2. The statement 
for arbitrary L then follows by generalization. 
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above threshold must vanish at the origin in order to satisfy orthogonality with respect to all 
such bound-state wave functions localized near the origin. In all of these cases the function 
bi{r) remains finite as r — > for any L. We conclude that for an effective field theory with 
contact interactions and attractive singular power-law interactions, we can take the cutoff 
momentum arbitrarily large for any L without producing a divergence in the coefficient 2 
of the short-range K matrix. 



VI. QUANTUM DEFECT THEORY AND THE MODIFIED EFFECTIVE 

RANGE EXPANSION 

Up to now we have been discussing the short-range phase shift of K matrix for scattering 
relative to the van der Waals wave functions Fl and Gl- For power-law interactions l/r a 
with a > 2, we also have the option to define phase shifts relative to the Bessel functions of 
the free wave equation. The problem though is that the usual effective range expansion, 



k 



2L- 



A cot 8 L , d {p) = h \r L k 2 

ClL 2 



D 

n=0 



(6.1) 



is spoiled by nonanalytic terms as a function of k 2 . For the van der Waals interactions in 
the L = channel, the leading nonanalytic term is proportional to fc 3 , and so the scattering 
parameters ao and tq are well defined, but the shape parameters V\/ are not. For L = 1 
the leading nonanalytic term is proportional to k and so only the scattering length a\ is well 
defined. For L > 2 none of the low-energy scattering parameters are well defined. To resolve 
these problems, a modified form of the effective range expansion is used which is known as 
quantum-defect theory [4"oT - T4"?] . 

In quantum-defect theory for attractive 1/r 6 potentials, one defines an offset for the phase 
shift [22], 

riL=\{v-v Q ). (6.2) 
The modified effective range expansion is then 



k 2L+1 cot (5 L + 2 VL ) 



-^ + \Rl k 2 + 0(kHnk) 



(6.3) 



where Al and Rl are the generalized scattering length and effective range parameters. These 
definitions coincide with the usual scattering length ai for L = 0, 1 and the usual effective 
range Tl for L = 0. The generalized scattering length and effective range can be written in 
terms of the short-range K matrix parameters as 



2/Q2L+1 



A, 



2«+i[r(f + i)r(L + §)p 



K 



L,0 



(6.4) 



and 



Rj 



2 4L+2 r(f + irr(L + j) ^, 

7r2(^L,0(-l) L -l) 



3\ 2 0-2I/-1 



AL 2 + AL - 3 



-K 



L.2 



(6.5) 



From these results we see that the short-range parameter Kl^ appears in combination with 
/3q. But in nearly all single- channel scatterings between pairs of alkali- metal atoms, from 
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the following equation, 

Pi 



K 



L.2 



1 + 



2^r(f + i) 2 r(L + |) 



3\ 2 



4L 2 + 4L - 3 

2^r(f + i)^r(L + l) 



A, 



-I -2 



vr 2 /3 2L+1 



2 4, + i r(f + i rr(L + § y 



2«2L+1 



7T 2 /?, 



(6.6) 



one quantitatively finds that Kl,2 is at least one order of magnitude smaller than /3|. This 
separation of scales is the reason for the approximate universality found in Refs. 



The dominance of /3| n over the sub leading coefficients K^n i n Eq. (4.1) for n > 1 holds 



for nearly all cases of single-channel scattering between alkali-metal atoms [2IJ HB]. This 
phenomenological fact explains the absence of short-distance length scales in the universality 



relation in Eq. (1.2). Furthermore, Gao has shown that when short-range interactions arise 
from a repulsive central potential, the fact that the K matrix is nearly independent of 



energy means that the K matrix is also nearly independent of angular momentum L 
This produces a surprisingly rich class of universal physics for single- channel van der Waals 
interactions where K^^n is negligible compared to for all L, and Kl Q is approximately 
the same for all L. Therefore /3q and the s-wave scattering length will determine, to a good 
approximation, the threshold scattering behavior for all values of L. 



VII. CAUSAL RANGE FOR SINGLE-CHANNEL SCATTERING 

We have shown that for negative < 0, the range R of the short-range interaction 
can be taken all the way down to zero. But when Kl i is positive, there is a constraint on 



R and we use Eq. (4.18) to determine a minimum value for R. We call this minimum range 



the causal range R", and we determine R b as the solution to the equation 

K L , 2 = b L (R b ). (7.1) 
As pointed out in Ref. [33J, one can show a priori that bi(r) is a monotonically increasing 



function of r. Therefore, if a real solution to Eq. (7.1) exists, then it is unique. If, however, 
there is no real solution, then there is no constraint on the interaction range and we define 
R h to be zero. For an effective field theory with contact interactions and van der Waals tail, 
the cutoff momentum can be made as large as ~ 1/R b before the causality bound is violated. 

In the following analysis we extract the single-channel s-wave effective range parameters 
ao and ro for several different pairs of alkali-metal atoms 7 Li, 23 Na, and 133 Cs in singlet and 
triplet channels. The data is taken from Refs. [201 12TI I49H5T] . The reduced masses for 7 Li 2 , 
23 Na 2 , and 133 Cs 2 are \i = 6394.7,20954,121100 (a.u.), respectively. The van der Waals 
coupling constants for 7 Li 2 , 23 Na 2 , and 133 Cs 2 are C 6 = 1388, 1472, 7020 (a.u.). We calculate 



the corresponding i^-matrix parameters using Eq. (6.4) and Eq. (6.6) and then compute the 
resulting causal ranges. We recall that for L = we simply have Ao = ao and Rq = ro. 
The results for the scattering parameters and causal ranges are given in columns II, V and 
VI of Table [1} The discrepancies in R are due to the fact that in the analytic studies in 
Refs. [201 EI] Ko,2 is neglected, while the numerical calculations of Refs. [49T15T] include the 
short-range contribution from i^o,2- 
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TABLE I: Scattering parameters and causal ranges for s-wave scattering of 7 Li, 23 Na, and 133 Cs pairs. 
The scattering data collection is taken from Ref. |20J. In columns I and IV the scattering data for 7 Li are 
from Ref. [50], the scattering data for 23 Na are from Refs. [49J [50] , and data for 133 Cs are from Ref. [51] . 

In column III the effective range parameters, Rq, are calculated analytically in Refs. [20 ] 121 ) . In column 
IV, the Rq are obtained from numerical calculations. The scattering parameters in columns II and V are 

calculated using Eqs. (6.4) and (6.6), and the causal ranges in column VI are obtained from Eq. (7.1). 







I 


II 


III 


IV 


V 


VI 


Atoms 


State 




A 


#0,0 


i?0 


R 


#0,2 


R b 


7 Li- 7 Li 


^9 


64.9097 


36.9 


-5.282 


66.3 


66.5 


2 - 124 


7 


- 19 


7 Li- 7 Li 


3y 


64.9097 


-17.2 


0.643 


1006.3 


1014.8 


0-17 


3 


- 25 


23 Na- 23 Na 


^9 


88.624 


34.936 


5.705 


187.317 


187.5 


0-86 


4 


- 20 


23 Na- 23 Na 


3y 


88.624 


77.286 


-1.213 


62.3756 


62.5 


2-13 


16 


- 24 


133 Cs _133 Cs 


^9 


203.62 


68.216 


3.365 


624.013 


624.55 


- 146 


7 


- 45 



In column V of Table [TJ we present an approximate range for K 2 for each atomic pair 
using the values for R in columns III and IV. Since K ^ is positive, we cannot go all the 
way to the zero-range limit. However, in each case Kq^ is at least one order of magnitude 
smaller than /3|. 2 Although we cannot take the zero-range limit, the causal ranges are small 
in comparison to /3q. In each case R b is less than one-third the size of (3$. Hence one can 
probe these interactions in a van der Waals effective field theory with cutoff momentum up 
to roughly three times l//3 6 without violating the causality bound. 



VIII. CAUSAL RANGE NEAR A MAGNETIC FESHBACH RESONANCE 



In Ref. [52] the multichannel problem of scattering around a magnetic Feshbach resonance 
is reduced to a description by an effective single- channel K matrix that depends on the 
applied magnetic field B. The behavior around the resonance is described by several 
parameters. B 0t i is the position of the resonance, while g res parametrizes the width of the 
Feshbach resonance. K^ g is a background value for the K matrix, and the scale cIb,l is 
introduced to define a dimensionless magnetic field. We write the effective single-channel K 
matrix as 



Kf(k, B) = -Kl 



with 



bg 



1 + 



9t 



k 2 Pl-g TeB (B s + 1)\ 



(B-B 



0,L) 



d 



B,L 



•1) 



(8.2) 



The parametrization given above corresponds to Eq. (18) in Ref. [52] ■ Note that we have 
changed the notation slightly and are using a diff erent sign convention. 

By expanding the right-hand side of Eq. (8.1) in k 2 , it is straightforward to determine 



2 Note that Kl2 h as the dimension of an area (in the appropriate atomic units). 
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the K matrix expansion parameters Klo and Kl,2- A short calculation yields that 



K 



eff 
L,0 



1 + 



1 



K eS _ 



g ieB (B s + r 



As noted in Ref. [52], the parameters K^ s and g res are constrained by the condition 



A^ rcs <0 



1.3) 



(8.4) 



(8.5) 



From this we directly see that K^ 2 given by Eq. (8.4) is always negative. From the causality 



bound in Eq. (4.18) it follows that where this effective single- channel description is applicable 



and correctly captures the entire energy dependence of the short-range K matrix, the causal 
range will be zero when the interaction is tuned close to a Feshbach resonance. 



IX. SUMMARY AND DISCUSSION 

In this paper we have analyzed two-body scattering with arbitrary short-range interac- 
tions plus an attractive 1/r 6 tail. We derived the constraints of causality and unitarity for 
the short-range K matrix, 



(9.1) 



n=0 



For any r larger than the range of the short-range interactions, R, we find that Kl i2 satisfies 
the upper bound 

K L , 2 < b L (r) , (9.2) 



where bAr) is 



b L (r) = ^W[f L>2 , f L>0 )(r) + lKl W[g L>2 , g Lj0 ](r) 



TV 



K L , W[g L>2 , f L , ](r) + W[f L>2 ,g Lfi ](r) , 



(9.3) 



and the Wronksians are given in Eq. (4.15), (4.16), and (4.17). 



In clear contrast with the case for finite-range interactions only [381 139] or with Coulomb 
tails [33], the function bi{r) does not diverge but rather vanishes as r — > for all L. When 
Kl,2 ^ 0, there is no constraint derived from causality and unitarity that prevents the use 
of an effective field theory with zero-range contact interactions plus an attractive 1/r 6 tail. 
This holds true for any angular momentum value L. For the phenomenologically important 
case of a multichannel system near a magnetic Feshbach resonance, the effective value for 
Kl j2 is negative and so the short-range interaction can be taken to have zero range. 

The van der Waals interaction is qualitatively different from the Coulomb interaction 



where 69' 



oulomb / 



diverges for L > 1. The key difference is that both van der Waals wave 



functions Fi and Gl vanish at the origin. This phenomenon also occurs for an attractive l/r a 
interaction for a > 2 in any number of spatial dimensions. It is also valid for an attractive 
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1/r 2 interaction when the coupling constant is strong enough to form bound states. For an 
effective field theory with contact interactions and attractive singular power-law tail, the 
cutoff momentum can be made arbitrarily large for any L without producing a divergence 
in the coefficient of the short-range K matrix. 

When 2 is positive, there is a lower bound on the range of the short-range interactions. 
We define the causal range R b as this minimum value for the range, given by the condition 

K L , 2 = b L (R b ) . (9.4) 

We have analyzed several examples of s-wave scattering in alkali-metal atoms in Table |T} 
We find that the Kl,2 is at least one order of magnitude smaller than /3f . As a result we 
find that the causal ranges are small in comparison with /3q. 

In summary, we find that /3q dominates over distance scales parametrizing the short-range 
interactions. The origin of this van der Waals universality can be explained by two facts. 
The first fact is the phenomeno logical observation that, in single-channel scattering between 
alkali-metal atoms, there is a significant separation between the typical length scales of 
the short-distance physics and /3q. This can be seen by the small size of the short-range 
parameter 2 compared with As Gao has shown, this also leads to the approximate 
universal relation that Kl is the same for all L [48J. Therefore, to a good approximation, 
/3q and the s-wave scattering length will determine the threshold scattering behavior for all 
values of L. For the multichannel case near a magnetic Feshbach resonance, we find that 
the effective KJ^ 2 is no longer negligible. However, K^ 2 is negative, and this means that 
there is no constraint from causality preventing the zero-range limit for the short- distance 
interactions. 

The second fact underlying the van der Waals universality is that the zero-range limit 
of short-distance interactions is well behaved with regard to scattering near threshold. We 
note, however, that there is still no scale-invariant limit for L > 1 since the effective range 
parameter will diverge to negative infinity as (3§ goes to z ero. This can be seen from the 



(3q 2L+1 behavior with negative coefficient for L > 1 in Eq. (6.5). 

The analysis in this paper should be useful in developing an effective field theory with 
an attractive 1/r 6 tail and contact interactions. Similarly, one can also construct effective 
field theories for other attractive singular potentials l/r a for a > 2. These effective field 
theories could be used to investigate the Efimov effect and other low-energy phenomena in 
a model-independent way. 
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Appendix A: Van der Waals wave functions 



In this section we derive the van der Waals wave functions and Gl, following the steps 
in Ref. [33]. We first redefine the radial function as Ul(t) = y/r^Z(p s ). This rearrangement 



puts Eq. (3.9) into the form of an inhomogeneous Bessel equation, 

„ d 2 



Ps 



dp\ 



+ ps 



d 

dp, 



with 



2 i 2 
- ^0 + Ps 



-(2L + 1). 



Z(ps 



K Z{ Ps 
8 Ps 



(Al) 



The idea, introduced in Ref. [53], is now to consider Z v (p s ) as a series expansion of solutions, 



(A2) 



n=0 



and to use perturbation theory to obtain a solution for Z v (p s ). Substituting Eq. (A2) into 
Eq. (Al) leads to an infinite number of differential equations, 

C V0 ^\ Ps ) +k] 



op s 



0. 



(A3) 



The zeroth-order differential equation is homogenous, while all other orders are inhomoge- 
neous. This procedure generates a secular perturbation in all inhomogeneous differential 
equations as well as driving terms. The secular terms here refer to the solutions of the 
zeroth-order differential equation, which are Bessel functions. 

Following Ref. [33], we introduce a function Z u (p s ) which has an expansion in terms of 
Bessel functions with momentum-dependent coefficients, 



Z u {p s ) = b m (k s ) J u+m (p s ) . 



(A4) 



m=— oo 



We insert this as an ansatz into Eq. (Al) with v yet to be determined. Here J n denotes 
collectively the Bessel functions of the first and second kind, J n and N n . Substitution of 



Eq. (|A4|) into Eq. (|A1|) yields a three-term recurrence relation for the b m functions with 

hi 



—oo < m < oo, 



[(u + m) -is }b m (k s ) + 



ip/ —rzb m -i(k s ) + — f — — — b m+ i(k s 

16(m + v — 1) 16(m + v + 1) 



0. 



(A5) 



Solving these equations for b m {k s ) yields 



bm{ks 



•l) r 



2m 



T(u)T(u-uq + 1)T(u + u q + 1) 
T(v + m)Y{y — uq + m + l)T(u + vq + m + 1 



c m {u) (A6) 
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and 



' h 



2 m 



Y{v — m + l)T(z/ — vq — m)Y{y + uq — m) 
i> + l)i>-z/ )i> + ^) 



Cm. V) 



(A7) 



for m > 0. The functions c m (±z/) are defined as 



m— 1 



c m (±z/) = fj Q(±z/ + s) 6 (A; S ) 



(A8) 



s=0 



where Q(^) is given by 

Qiy) = — 



A' 



(A9) 



16(1/ + + l) 2 - z/ 2 ](z/ + 2)[(u + 2) 2 - v*] 

The coefficient bo(k s ) only determines the overall normalization and is simply set to one in 
the following. Equation (A5) for m = determines the shift v in the order of the Bessel 
functions. We determine v using the constraint 

[ 0) Wu(u-l)[(u-iy-ui} Ks 16V(z/+l)[(z/ + l)2-Z/2] Ks lA1UJ 



In general there are several roots which become complex beyond a critical scaled momentum 
k s , and one must be careful to choose the physical solution. For a detailed discussion of this 
point, see Refs. [221 Hi]. 

Choosing either J n = J n or J n = N n already yields a pair of linearly independent 
solutions. However, in order to get a pair with energy-independent normalization as r s — > 
(which ensures analyticity in the energy), we furthermore define 

oo oo 

X L (k s )=cos VL (-l) m M*.)-sin77£ ^ (-l) m & 2m+ i(^) (All) 



m=— oo 



m=— oo 



and 



with 



n(fc.) = sin»fc X] (-1)^(^+008^ J2 (-i) m hm + i(ks) , (A12) 



Combining everything, we arrive at the van der Waals wave functions, 



F L (k,r) 



G L (k,r) 



,1/2 



Xl(k s ) + Y£(k s ) 

1/2 

ry 



^l(^s) ^2 b m (k s ) J u+m (p s ) - Y L (k s ) 22 b m (k s ) N u+m (p s ) 

m=— oo m=— oo 

(A13) 

oo oo 

X L {k s ) b m (k s ) N u+m (ps) + Yl(£;) ^ b m (k s ) J v+m {p s ) 

(A14) 



m=— oo 



m=— oo 
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Appendix B: Low-energy expansions 



In this appendix we expand all functions relating to the van der Waals wave functions in 
powers of momentum. We first consider u, the shift in the order of the Bessel functions in 



(Bl) 



Eq. ( |330| ) and Eq. ( |3J1[ ). Using Eq. < \M0\ ) in Appendix |AJ we find 

3 



v = i/ - 



2V (4z, 2 -l)K 2 -l) 



K + o(K 



where v = (2L + l)/4. Using the expansion in Eqs. (A6), (A7), and (A9), we get 

2m 



:-iy 



r(z/ )r(2z/ + 1) 



and 



i\k s 



ml T(u + m)T(2vQ + m + 1) \ 4 
r(z/ — m + l)r(2i/ — m) f k 



+ 0(k 



2m+2\ 



m!r(i/ + l)r(2i/ ) 



2i» 



+ 0(A; 



2m+2\ 



for m > 0. Substituting these expressions into Eqs. (All) and (A12) we obtain 



Y L (k a 



X L (k s ) = l + 0(kt), 

r(u )T(2u - 1) , r(i/ )r(2i/ + i; 



+ 



i>o + i)r(2z/ ) r^ + i)r(2i/ + 2) 



(B2) 

(B3) 

(B4) 
(B5) 
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